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We study the coupling induced destabilization in an array of identical oscillators 
coupled in a ring structure where the number of oscillators in the ring is large. The 
coupling structure includes different types of interactions with several next neighbors. 
We derive an amplitude equation of Ginzburg-Landau type, which describes the 
destabilization of a uniform stationary state in a ring with a large number of nodes. 
Applying these results to unidirectionally coupled Duffing oscillators, we explain the 
phenomenon of a fast transition to chaos, which has been numerically observed in 
such systems. More specifically, the transition to chaos occurs on an interval of a 
generic control parameter that scales as the inverse square of the size of the ring, i.e. 
for sufficiently large system, we observe practically an immediate transition to chaos. 
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The understanding of the dynamical behavior of networks of coupled oscillators 
can contribute to the explanation of various collective phenomena that can be 
observed in coupled systems in biology, economy, or physics 1 3 . Discrete media, 
e.g. in neural systems, can exhibit complex coupling structures including feed- 
forward loops, large coupling ranges, and interaction mechanisms of different 
kind. In such systems, the influence of the coupling can transform the simple 
dynamics of a single oscillator into complicated spatio-temporal structures of the 
network dynamics. Amplitude equations of Ginzburg-Landau type have been 
a powerful tool for a universal description of spatio-temporal dynamics and 
pattern formation in continuous media. In this paper, we apply this technique 
to a large class of coupled oscillator systems with a ring structure, where we 
assume that the number of oscillators in the ring is large. In addition, we use 
the amplitude equation together with the corresponding scaling laws to explain 
the coupling induced fast transition from homogeneous stationary behavior to 
high dimensional chaos, which can be observed in certain coupled oscillator 
systems. 



I. INTRODUCTION 

Networks of coupled oscillators have been the subject of extensive research in the last 
decade. Coupled systems can display a huge variety of dynamical phenomena, starting from 
synchronization phenomena in various types of inhomogeneous or irregular networks, up to 
complex collective behavior, such as for example various forms of phase transitions, traveling 
waves^-, phase- locked patterns, amplitude death states 7 , or so called chimera states that 
display a regular pattern of coherent and incoherent motion^"—. Of particular interest are 
situations, where complex spatio-temporal structures can emerge in regular arrays of identi- 
cal units induced only by the coupling interaction. In many cases, the resulting phenomena 
differ substantially from corresponding situations in continuous medial and depend strongly 
on the underlying network topology. 

Our specific interest is in the emergence of spatio-temporal structures in a homogeneous 
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array of identical units that have a stable uniform equilibrium at which the coupling vanishes. 
As a classical paradigm, the Turing instability gives an example of a coupling induced 
instability in such a setting. This phenomenon has of course a direct counterpart in the 
discrete setting, but it turns out that there appear some genuinely new phenomena. In-i2r— 
it has been shown that also in a ring of unidirectionally coupled oscillators, i.e. in a purely 
convective setting, the Eckhaus scenario of coexisting diffusive patterns can be observed. 
ImM it has been shown that Duffing oscillators coupled in the same way, exhibit a complex 
transition to spatio-temporal chaos. In this paper we develop a general theoretical framework 
for such phenomena in large arrays. We derive an amplitude equation of Ginzburg-Landau 
type that governs the local dynamics close to the destabilization of the uniform steady state. 
It resembles several features that are already well known in the context of reaction-diffusion 
systems^ 1 ^. But in contrast to these results, it can be applied to much more general 
coupling mechanisms, including also the case of unidirectional and ant i- diffusive interaction 
and allows also for a mixture of such interactions with several next neighbors. Such an 
interplay of attractive and repulsive coupling is an essential feature of neural systems. As 
a specific feature, the convective part will appear in the amplitude equation as a rotation 
of the coordinates in an intermediate time scale that is faster than the diffusive processes 
described by the Ginzburg-Landau equation. 

Having deduced the amplitude equation and the corresponding scaling laws in terms of 
the number of oscillators, which is assumed to tend to infinity, we use this theory for the 
explanation of a specific phenomenon that has been reported in^ for a ring of unidirectionally 
coupled Duffing oscillators. There, it has been shown numerically that for a large number of 
oscillators N, there is an almost immediate transition from homogeneous stationary behavior 
to high-dimensional chaos. Based on our amplitude equation, we argue that in such systems, 
one can expect generically that such a transition occurs within a parameter interval of the 
size 1/N 2 . We consider a generic case, where the control parameter enters already the linear 
parts of the dynamical equations, e.g. a diffusive coupling strength. Finally, we demonstrate 
this phenomenon by a numerical example, where we also evaluate the scaling behavior of 
the parameter interval where the transition to chaos takes place for an increasing number 
of oscillators. 



3 



II. MODEL EQUATION, SPECTRAL CONDITIONS, AND NOTATIONS 
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Figure II. 1. An example of a ring of N coupled oscillators. Apart from the self-coupling Mq, each 
oscillator yj is also coupled with yj+2 (M2) as well as yj-i (M_i). See Eq. (jll.ip for the equation 
of motion. 

We are interested in a system of N identical coupled oscillators that has an uniform 
equilibrium, where the coupling vanishes. The coupling network is organized in a ring 
structure, where interactions of several next neighbors are possible. Such systems can be 
written in general form as 

Vi = ^2 M m{p)yj+m + h(yj,y j+1 , . ..y j+N ;p), (II.l) 

m 

where yj G M", j = 1, . . . , N describes the state of the j-th oscillator. The ring structure 
is induced by periodic boundary conditions, i.e. all indexes have to be considered modulo 
N. The linear part of the dynamics is given by the n x n matrices M m (p), m = 1, . . . , N, 
depending on the bifurcation parameter p, which account for the coupling to the m-th 
neighbor; in particular M Q (p) describes the linear part of the local dynamics (self-coupling). 
The nonlinear part h, again including a local dependence and a dependence on the m-th 
neighbor, should vanish at the origin h(0, . . . , 0;p) = and have also zero derivatives there. 
Note that this system is symmetric (equivariant) with respect to index shift. Figure III. II 
illustrates an example with self coupling and coupling to the neighbor on the left and to 
the second neighbor on the right. The specific form of (III. II) also implies that the coupling 



vanishes at the equilibrium y\ — ■ ■ ■ = y^ = 0, which is true e.g. when the coupling is a 
function of the difference yj — y m for any two coupled oscillators j, m. 



The characteristic equation for the linearization at the zero equilibrium of (1II.1|) can be 
factorized as 



X (p, A, e^/ N ) 



det 



Xld-J2e im2nj/N M m {p) 



0. 



where Id denotes the identity matrix in M n and the index j = 1,2, ... ,N accounts for 
the iV-th roots of unity that appear as the eigenvalues of the circular coupling structure^. 
Following the approach in^^, we replace for large iV the discrete numbers 2irj/N by a 
continuous parameter ip>, and obtain the asymptotic continuous spectrum 



A, 



A G C : X {P, A, e llp ) = det 



AId-^e^M m (p) 



0, (f G [0, 2tt) 



(II.2) 



which contains all eigenvalues and which for large iV is covered densely by the eigenvalues. 
Since the expression (III. 2ft is periodic in <p, the asymptotic continuous spectrum A p has 
generically the form of one or several closed curves \ p (ip) in the complex plane, parametrized 
by tp. 

At the bifurcation value p = 0, we assume that the asymptotic continuous spectrum 
touches the imaginary axis at at some point iuo (see Fig. III. 21) . i.e. the following conditions 
are fulfilled 



Ao(<£o) = iuo, 



(Vc 



IK X , 



Ki G 



(II.3) 



The first condition from (III. 3[) means that the point iuo belongs to the asymptotic continuous 
spectrum, while the second condition from (III. 3[) guarantees that the real part K(A(</?)) is 
tangent to the zero axis at (p = ip (see Fig. III.2j) . This tangency condition is quite natural 
and describes the condition for the destabilization (or bifurcation) of the zero solution a 
large ring of coupled oscillators. Indeed, if the spectrum A p is contained in the left half 
of the complex plane with Re A < 0, then the uniform equilibrium y\ — ■ ■ ■ — — is 
asymptotically stable. As soon as A p crosses the imaginary axis, it becomes unstable for 
sufficiently large N. 

Before we present our main result, we now introduce some useful notations and formulate 
a technical lemma, which follows from the bifurcation conditions ( 1II.3|) . With Vq and V\ we 
denote the eigenvector and the adjoint eigenvector to the critical eigenvalue Ao(y?o) = i^o, 
which we assumed in (|II.3|) to exist for p = 0. Moreover, it will be convenient to denote by 



£o = $>"^°M m (0), 




Figure II. 2. Asymptotic continuous spectrum A p at the destabilization; schematically. 

L 1 = ^me im ^M m (0), 

m 
m 

the "moments" of the coupling matrices M m . In this notation, the equations for the eigen- 
vectors v and V\ read as 

[iwold - L ] v = 0, (II.4) 



[-iw Id - L* ] v x = 0. (II.5) 

These vectors can be normalized as 

\v \ 2 = 1, (« ,vi) = l. (II.6) 

Finally, we expand the matrices M m (p) with respect to the parameter p and write them 
them as M m (0) + pif m + O (p 2 ). In this way, we can define 



L K = J2e imw K m . 



Lemma 1. Assume that (po is a regular point of the asymptotic continuous spectrum lill.ty) . 



such that Ao(y?) exists and is locally differentiable in a small neighborhood of ip . Further, 
let the bifurcation condition UI. 3\) hold. Then 

(Lif ,fi) = Ki, (11.7) 

The proof of this Lemma will be given together with the proof of the main result in 
Section HVl 
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III. MAIN RESULT: REDUCTION TO GINZBURG-LANDAU EQUATION 



In this section, we present an amplitude equation that describes for system (IH.ip the 
dynamics close to the destabilization threshold in the limit of large N. We perform a 
limiting procedure, where we assume that a fixed number R of neighbors is involved in the 
coupling while the total number of oscillators N tends to infinity. Also the nonlinearity h 
is assumed to depend on R next neighbors only Hence, in the sequel all summation will 
extend over m G {—R, • • • , R} and the nonlinearity will be written as h(yj-R, . . . , yj+R,p), 
independently on N. In this way, the coupling becomes local in the limit N — > oo, and 
coupling terms will be approximated by derivatives of the amplitude. Consequently, our 
results will be valid for large rings, where the coupling range is small compared to the total 
size. It will be shown that the amplitude equation has the form of a complex Ginzburg- 
Landau equation with periodic boundary conditions. The derivation is quite technical and 
we give here in Proposition 2 only the main assertions. The proof is deferred to Section IIVI 

Proposition 2. Assume that the bifurcation condition ill. 3) holds and tpo is a regular point 
of the asymptotic continuous spectrum A . Additionally, let the points ±3iu)o not belong 
to the asymptotic continuous spectrum with tp — (fo (nonresonance condition). Let the 
nonlinearity h be of third order. Introduce the small parameter 

1 

£= N 

and apply the multiple scale ansatz 

yj (t) = ee^ t+i ™v A{T l: x u T 2 ) + e 3 e 3i( - u} ° t+ ^v 2 A 3 (T 1 ,x 1 , T 2 ) + c.c, (III.l) 

with the amplitude A G C depending on the rescaled coordinates 7\ = et, T 2 = e 2 t, and 
x\ = ej (c.c. denotes complex conjugated terms, v 2 G C n ) to the following system 

R 

y 3 = {M m (0)+e 2 rK m + O(s i ))y J+m + h(y J . R ,...y j+R ;p), (III.2) 

m=—R 

with the rescaled parameter p = e 2 r and j = 1, . . . , N with periodic boundary conditions. 
Then, the solvability conditions up to the order e 3 imply the following partial differential 
equation of Ginzburg- Landau type 

dr 2 u = rn 2 u + -^-d^u + (u\u\ 2 (III. 3) 
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with periodic boundary conditions 



«(e,T 2 ) = u(e+i,r : 



2), 



where u(£, T 2 ) with £ G [0, 1] zs related to the amplitude A by 

A{T u x u T 2 ) = u(k 1 T 1 + x 1 ,T 2 ). 



(III.4) 



The coefficient K\is given by ill. 3\) or UI. 7j) . and 



n 2 = (L K v , Vx) , k 3 = (L 2 v , Vx) . 

Finally, the coefficient of the nonlinearity ( and the vector v 2 G C n have to be determined 
by the nonlinearity h according to UV. 1\) and UV.ty . 



According to this proposition, small solutions of a coupled system of the form 2 j) that 
has a parameter close to criticality, i.e. p = 0(e 2 ) can be approximated in the form (1III.1I) . 
where the amplitude A is related to a solution of the Ginzburg-Landau equation 3 j) via 
flIII.4p . Note that the relation (1III.4P introduces a rotating frame on the timescale T\ = et 
with rotation velocity K\. The time evolution of the Ginzburg-Landau equation enters only 
on the slowest time scale T 2 = e 2 t. 

Remark 3. For the case of a symmetric coupling (M& = M_fc) the amplitude equation (1111.30 
has all real coefficients Ki, k 2 and k 3 and the rotation velocity vanishes Kx = 0. 



IV. PROOFS 



Proof of Lemma 1. We have to show that Kx, defined by (111.30 . can be calculated as given 
in (III.7p . To this end, we differentiate the eigenvalue equation (cf. (111.4j) ) 



[A (<p) ld-L ]v((p) = 0, 



with respect to ip: 







Id— AoM - i me im ^M m (0) 



v(<p) + 



AoMld-^e^M m (0) 



Evaluating the obtained expression at ip = <po, \(<po) = iu>o, and v(ipo) = Vq, we obtain 

dv 

[IdiKx - iLx] v + [iu ld - L ] q^(Vo) = 0. 
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The projection onto Vq gives 

/ dv 
«i (v ,vi) - i {Liv Q ,vx) + (Hu ld - L ] — (</?o),^i ) = 0, 

which implies 

/ dv 



Taking into account (111.51) . we obtain the relation ( 1II.7[) . which proves the Lemma. 



Proof of Proposition 2. Substituting the multiple scale ansatz flHI.ll) into (IIII.2I) . we 

obtain 

at v ' 
= eJ2 ( M m(°) + £ 2 rK m ) e iuiot+ivoj e imipo VoA(T u x x + me, T 2 ) 

m 

+e 3 J2 ( M -(0) + e 2 rK m ) e^+^V^ ^ 3 ^, ^ + ^ ^ + c c 

rn 

+h(yj-R, • • • , e 2 r) 

Dividing the obtained equation by e, expanding necessary arguments of A, we obtain up to 
the terms of the order e 2 (the complex conjugated terms are omitted here for brevity) 

e^+^vo (itu A + ed Tl A + e 2 d T2 A) + Ziuj e 2 e 3l{uJot+wj) v 2 A 3 
= e iuot+l * oj ( M ™(°) + £ 2 rK m ) e mvo v (a + med Xl A + ~m 2 e 2 d 2 Xl A) 

+e 2 M m {0)e 3i{uiot+ ^ j) e 3imvo v 2 A 3 

m 

+£ 2 e iu t+iwi A |^|2 ^ ^ + £ 2 e 3i(uot+<p j) A 3 h2 ^ 

where hi (vq) and h 2 (vq) are determined by the leading terms in the expansion of the non- 
linearity. Note that due to our assumption, h is of third order and the leading order terms 
are given by expanding h(y.j-R, . . . , Uj+R, 0) in the homogeneous state y m = av^ + cc, m = 
j — R, . . . , j + R with respect to a £ C The solvability condition requires that different 
harmonics as well as different orders of e up to e 2 are equal. We start with the first harmonic. 
By equating the terms containing e lw o*+«w we obtain the following equation 

v (iu A + ed Tl A + e 2 d T2 A) = 
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= ( M ™(°) + e%rK m) e imtfo Vo U + med Xl A + Un 2 e 2 d 2 Xi A\ + e 2 A \A\ 2 h x (v ) 
Since it should be satisfied for all e, we first consider the e° equation 

iu Av = AL v , 

which holds according to the spectral condition flll.4p . The e 1 terms result into 



v d Tl A - LiV d Xl A = 0. 
Multiplication with vj and using ( III. TP from Lemma 1, we obtain 



9tiA — K\d xl A — 0. 
This will be accounted for by introducing the new amplitude u by 

tt(f , T 2 ) = u(k\T\ + x u T 2 ) = A(T 1 ,x 1 ,T 2 ) 
in a correspondingly rotating coordinate £ = KiTi + Xi. Finally, the e 2 terms result into 

v d T2 A = rAL K v + ^d 2 Xi AL 2 v + A \A\ 2 h x (v ) 

Note that the dependence on T\ does not show up in this equation. Hence, after multiplica- 
tion with v±, we can write it in terms of u as 



^3p2„. , /*-_ . I..I2 

2 



8t 2 U = TK 2 U + — d^u + I'" 
where 

C=(MVMi>- (iv.i) 

Finally, it is simple to check that the solvability of the terms for the third harmonic leads 
to the expression 



V-2 



3M)-^M m (0)e 



-1 



h 2 (v ) (IV.2) 



Here, the existence of a nonzero solution v 2 is guaranteed by the nonresonance condition. 
Indeed, since the points ±3iuq do not belong to the asymptotic continuous spectrum with 
tp = ifo, the matrix [3iuo — Yl m M m (0)e 3vmipo } is non-singular. Finally, note that the set of 
equations should be complemented by periodic boundary conditions in £, taking into account 
the ring structure of system flIII.2p . The proposition is proved. 
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V. APPLICATION: SCALING OF THE BIFURCATION TRANSITION 



It is well known that in the complex Ginzburg-Landau equation f]III.3j) a destabilization 
transition is possible, where in a series of bifurcations a homogeneous stationary state loses 
its stability and, via periodic and quasiperiodic motion, a regime of spatio-temporal chaos 
appears^"— . We will now use the result about the amplitude equation in order to trace 
back this destabilization transition to the original coupled oscillator system of the form 
(IH.ip . Specific attention will be paid to the scaling behavior of the parameter region where 
this transition takes place. 

Starting with the scale free amplitude equation 

0t 2 u = tk 2 u + -^-d^u + (u\u\ 2 (V.l) 

with periodic boundary conditions, we assume that at r = r a destabilization of the ho- 
mogeneous stationary state takes place and for further increasing parameter r a transition 
to spatio-temporal chaos can be observed. This transition takes place within some bounded 
interval Ar. Considering now the corresponding behavior of an oscillator system that is de- 
scribed by this amplitude equation, we have to introduce the scalings given in Proposition 
[2j In particular, according to the parameter scaling p = re 2 , the same transition takes place 
in the coupled oscillator system within a parameter interval 

Ap(N) = Are 2 = Ar^, 

of the corresponding parameter p of the oscillator system. Therefore, one can expect that 
a generic transition to chaos (or hyperchaos) in the system of N coupled oscillators occurs 
within a parameter interval that scales as 1/N 2 . In the following section, we present an 
example where the assumptions that we made above are satisfied and, consequently, in the 
oscillator system an almost immediate transition from stationary behavior to chaos can be 
observed. 



VI. EXAMPLE: A RING OF UNIDIRECTIONALLY COUPLED DUFFING 
OSCILLATORS 

In the previous section we discussed the scaling law for the transition to chaos in a ring of 
coupled oscillators. In this section we provide a numerical example to illustrate this scaling. 
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We use the autonomous Duffing oscillator described by the following second order ordinary 
differential equation 



y + dy + ay + y 3 = 0, (VI.l) 

where d and a are positive constants. System (IVI.ip is a single-well Duffing oscillator, which 
has a single equilibrium point at y — y — 0. Due to the presence of damping (d > 0) 
this equilibrium is an attractor for all initial conditions. We consider now a ring of TV 
such oscillators with a linear unidirectional coupling to the next neighbor. Introducing new 
coordinates x = y, z = y and the coupling into Eq. (IVI.ip . the equations of motion have 
the form 

X] = Zj ' (VI. 2) 

z j — clx j x a I k (yOb j *^ j ) ? 

where k is the coupling coefficient and indices are considered modulo N. A detailed anal- 
ysis of the behavior of this system was already presented in^. For at least three coupled 
oscillators and increasing coupling strength k, one can observe rich dynamics starting from 
periodic oscillations to hyperchaos. Here we focus our attention on the transition to chaos 
(hyperchaos) and its dependence on the number of oscillators. In our numerics, we used 
the fixed parameter values a = 0.1 and d = 0.3. We calculated the maximum Lyapunov 
exponents for increasing N and varying coupling coefficient k. based on this, we determined 
two values: kjj where the uniform stationary state loses its stability in a Hopf bifurcation, 
and kch where the transition to chaos takes place. A reliable computation of the Lyapunov 
exponents was possible only for moderate values of N {N < 30). For larger systems, we 
monitored the behavior in an appropriately chosen Poincare section in order to determine 
kn and kch- Fig. IVI.l( a) shows that for increasing N the distance between kn and kch 
decreases. Finally, for a large enough number of oscillators, the transition to chaos appears 
practically immediately after the Hopf bifurcation. In order to validate our conjecture about 
the 1/N 2 scaling, we plot in Fig. IVI.l( b) the scaled transition intervals 

k Re = (k Ch -k H )N 2 . 

The numerical results indicate that the scaled transition interval neither tends to zero nor 
diverges to infinity, which supports the scaling assumption. 
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Figure VI. 1. (a) Values of the coupling parameter kn at Hopf bifurcation (crosses) and kch a t the 
transition to chaos (circles) for increasing number of oscillators N. (b) Rescaled transition interval 
kR e (circles) versus N. 



In Fig. IVI.2I we present the Lyapunov spectrum for thirty coupled oscillators. Panel (a) 
shows that at kch = 0.164 two Lyapunov exponents become positive practically simultane- 
ously. A similar behavior has been observed for other large values of N. This shows that 
the transition to chaos and to hyperchaos occursin a large ring almost simultaneously. With 
further increasing coupling coefficient k (Fig. IVI.2f b)) one can observe that more and more 
Lyapunov exponents become positive which leads to a high dimensional chaos^ 1 ^. 




Figure VI. 2. Lyapunov exponents for N = 30 unidirectionally coupled Duffing oscillators (a) five 
largest Lyapunov exponents, ku and kch indicate the Hopf bifurcation and the transition to chaos, 
k € (0.139,0.17), (b) twenty largest Lyapunov exponents for k € (0.0,0.8). At k = 0.8, there are 
fourteen positive Lyapunov exponents. 
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VII. CONCLUSIONS 



The results of this paper are twofold. The first one is more theoretical: we have derived 
an amplitude equation for small amplitude solutions of a ring system of coupled oscillators. 
This equation has the form of a complex Ginzburg-Landau equation with periodic bound- 
ary conditions. Such an amplitude equation is typical for reaction-diffusion systems^ 1 ^ in 
continuous media and can also be found for systems with large delay 25 . Obviously, a sim- 
ilar behavior should be expected in discrete systems with diffusive coupling. However, our 
results show that in coupled oscillator systems the Ginzburg-Landau equation can also be 
used to describe the dynamics in systems with unidirectional, i.e. purely convective coupling 
or even with anti-diffusive interaction. In this sense, the class of coupled oscillators that we 
treated in this paper differs substantially from discrete analogs of the classical results for 
continuous media. 

For the second result, we applied the amplitude equation to the specific scenario of 
the transition to spatio-temporal chaos in such systems. In this way, we provided a general 
framework for the coupling induced transition to high dimensional chaos in coupled oscillator 
systems, that has been described for a specific example inM. In particular, we were able 
to show that the observed fast transition to chaos is a generic feature and follows a 1/N 2 
scaling law for systems with a large number N of oscillators. As a result of this scaling, 
one observes for large systems a practically immediate transition from a uniform stationary 
state to chaos. We illustrated this behavior by a numerical example, where the emergence 
of high dimensional chaos is demonstrated by corresponding Lyapunov spectra. 
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